Induced-charge electrophoresis near an insulating wall 



O 

o 
o 

Q 



C/3 



X3 
O 

o 



> 
m 

o 

O 



X 



Mustafa Sabri Kilic^ and Martin Z. Bazant^'^ 
Department of Mathematics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA and 
^ Physico-Chimie Theorique, Gulliver- CNRS, ESP CI, 10 rue Vauquelin, Paris 75005, France 

(Dated: February 2, 2008) 

Induced-charge electrophoresis (ICEP) has mostly been analyzed for asymmetric particles in an 
infinite fluid, but channel walls in real systems further break symmetry and lead to dielectrophoresis 
(DEP) in local field gradients. Zhao and Ban (Langmuir, 23, 2007, pp 4053) recently predicted that 
a metal (ideally polarizable) cylinder is repelled from an insulating wall in a DC field. We revisit this 
problem with an AC field and show that attraction to the wall sets in at high frequency and leads to 
an equilibrium distance, where DEP balances ICEP, although, in three dimensions, a metal sphere 
is repelled from the wall at all frequencies. This conclusion, however, does not apply to asymmetric 
particles. Consistent with the recent experiments of Gangwal et al. (|arXiv:0708.2417[l . we show that 
a metal/insulator Janus particle is always attracted to the wall in an AC field. The Janus particle 
tends to move toward its insulating end, perpendicular to the field, but ICEP torque rotates this end 
toward the wall. Under some conditions, the theory predicts steady translation along the wall with 
an equilibrium tilt angle, as seen in experiments, although more detailed modeling of the contact 
region of overlapping double layers is required. 
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I. INTRODUCTION 

Most theoretical work on electrophoresis has focused on 
spherical particles moving in an infinite fluid in response 
to a uniform applied electric field [3, [O, [13, HI]. Of 
course, experiments always involve finite geometries, and 
in some cases walls play a crucial role in electrophore- 
sis. The linear electrophoretic motion of symmetric (spher- 
ical or cylindrical) particles near insulating or dielectric 
walls H, [T^ 



or cylindrical)^ particles near ii 

20, ^2^1 and m bounded cavities or chan- 
nels |a[alllllii,i5., 21, 22, 24, 37, 48] has been analyzed ex- 
tensively. Depending on the geometry and the double-layer 
thickness, walls can either reduce or enhance the transla- 
tional velocity, and the rotational velocity can be opposite 
to the rolling typical of sedimention near a wall. The classi- 
cal analysis for thin double layers assumes "force-free" mo- 
tion driven by electro-osmotic slip alone, but recent work has 
shown that electrostatic forces can also be important near 
walls [igI . [45| . Heterogeneous particles with non-uniform 
shape and/or zeta potential exhibit more complicated bulk 



motion lU9|, l25l |26[ , which can also affect boundary inter- 
actions 42], especially if the particles are deformable, as in 
the case of chain-like biological molecules [3l| . 

In this article, we focus on the effect of nonlinear induced- 
charge electro-osmotic (ICEO) flows at polarizable surfaces, 
which are findirig many new applications in microfluidics and 
colloids d, [H, [4l[ . The canonical example of quadrupolar 
ICEO flow around a polarizable particle, first described by 
Murtsovkin [l^, [13] , involves fluid drawn in along field axis 
and expelled radially in the equatorial plane in an AC or 
DC field. Broken symmetries in this problem can gener- 
ally lead to hydrodynamic forces and motion induced-charge 
electrophoresis (ICEP), as well as electrical forces and mo- 




(a) Insulating wall 




(b) Conducting wall 



FIG. 1: Hydrodynamic forces on polarizable particles near (a) 
insulating and (b) unscreened conducting walls due to ICEO flows 



tion by dielectrophoresis (DEP). Such phenomena have only 
been analyzed for isolated asymmetric particles in an infi- 
nite fluid II, [40l . [4^ or in a dilute solution far from the 
walls [s^-ISo). In contrast, experiments demonstrating trans- 
lational ICEP motion have involved strong interactions with 
walls 11, 29], which remain to be explained. 

As shown in Figure [U it is easy to see that the quadrupo- 
lar ICEO flow around a polarizable particle typically causes 
attraction to unscreened conducting walls (perpendicular to 
the field) and repulsion from insulating walls (parallel to the 
field). The former effect of ICEP attraction to conduct- 
ing walls has not yet been analyzed; it may play a role in 
colloidal self assembly on electrodes applying AC voltages 
[3^ . [33I [ssl . [43I . [46j . This phenomenon is mainly understood 
in terms of electrohydrodynamic flows (what we would term 
"ICEO") induced on the electrodes, not the particles (typi- 
cally latex spheres), but ICEP could be important for more 
polarizable particles. 

The latter effect of ICEP repulsion from insulating walls 
has recently been analyzed by Zhao and Bau [47j in the case 
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of a two-dimensional ideally polarizable cylinder in a DC 
field. However, this phenomenon has not yet been confirmed 
experimentally. On the contrary, Gangwal et al fl| have 
recently observed that metallo-dielectric Janus particles are 
attracted to a glass wall, while undergoing ICEP motion par- 
allel to the wall and perpendicular to an applied AC field. 
It is not clear that the existing theory of ICEP can explain 
this surprising behavior. 

The objective of this work is to analyze the motion of 
three-dimensional polarizable particles near insulating walls 
in AC fields. As summarized in section [III we employ the 
standard low-voltage model in the thin double-layer approx- 
imation, following many authors [1, IH, IH, HO, |3| , including 
Zhao and Bau [47|. In section lllll we first analyze ideally po- 
larizable cylinders and spheres near a non-polarizable wall, 
which only experience forces perpendicular to the wall. In 
section IIVI we then study spherical metal/insulator Janus 
particles, which are half ideally polarizable and half non- 
polarizable. Due to their broken symmetry, the Janus parti- 
cles also experience ICEP and DEP torques, which strongly 
affect their dynamics near the wall. 



II. MATHEMATICAL MODEL 

A. Low Voltage Theory 

In this paper, we will consider either a cylindrical or a 
spherical particle of radius a in a semi-infinite electrolyte 
bounded by a plane. The distance between the center of the 
particle and the plane is denoted by h. In the absence an ap- 
plied electric field, we assume that the particle and the wall 
surfaces are uncharged. In addition, we will assume the elec- 
trolyte has a low Reynolds number, and impose Stokes equa- 
tions. We will assume that the thin double layer approxima- 
tion holds and the bulk electrolyte remains electroneutral, 
which is the case when the Dcbyc length 



A, 



ekT 



is much smaller than the characteristic length scale (in our 
case, a). The Debye length is typically ranges between 1 — 
lOOnm, and colloidal particles are usually in the /i range, 
therefore thin double layer approximation holds for most of 
the time. 

Then the general equations consist of the Laplace's 



eV^(j) = 



and Stokes equations 



r/V^u 
V • u 



Vp 




(1) 



where (j) is the electrostatic potential and e the permittivity, 
T] the viscosity of the electrolyte, u the velocity field and 



p the pressure. The wall boundary z = is an insulator, 
satisfying 

n • V</> = 

whereas the particle surface, being polarizable, acts 
pacitor in the thin double layer limit 

with q being the surface charge density on the particle, k 
the conductivity of the bulk electrolyte. Far away from the 
particle, an electric field 



V0 - Eoo = Eaoi, |X| 



Eooit) may be 



is applied. In general, the magnitude 
time dependent. 

The Stokes equations are supplied by the no-slip condi- 
tions on the wall and the Smoluchowski's electrokinetic slip 
formula on the particle 

£ 

V 

where C is the potential difference between the surface and 
the bulk. Later, we will take into account the steric effects of 
ions and show how this formula needs to be modified. Lastly, 
we assume flow vanishes at the infinities. 

So far, the equations are complete except for a constitutive 
relation between C and q. The linear theory asserts that 

, = -fc 

Ad 

although the complete problem is still nonlinear in this case 
because of the quadratic slip formula ([S]). In this paper, we 
will study only the linear theory, but the calculations can be 
repeated with (sometimes more accurate) nonlinear theories 
like that of Poisson-Boltzmann. 



B. Force and Torque on the Particle 

Throughout this paper, we will assume that the particle is 
fixed and calculate the forces on the particle. For the case of 
a moving particle, the slip velocity needs to be modified to 
account for the motion of the particle surface. 

The total force and torque on any volume of the fluid are 
conveniently given in terms of the stress tensor, a, by 



F = / n ■ adA 



/ r X (n • a)dA 
Jan 



(2) 
(3) 



The stress tensor contains from electrical and viscous stresses 
on the fluid, a = a m + <^H, where 

aM = e[EE - ^EH] 



an = -pI + T] (Vu + (Vu)^ 



are the Maxwell and hydrodynamic stress tensors, respec- 
tively. 



C. Particle Dynamics 

In order to calculate the movement of a colloidal particle, 
we need to find a translational velocity U, and a rotational 
velocity rj such that the net force on the particle is zero, 
when the slip velocity is modified by taking into account the 
velocities U and ft. In other words, we are seeking U and fl 
such that the problem ([T|) with boundary condition 



-U + r X J7 



yields F = and T = 0. 

Since Stokes problem is linear, there is a linear relation- 
ship between the translational and rotational motion of the 
particle and the resulting force and torque exerted on it by 
the fiuid. Let us denote this relationship by 



M 



The viscous hydrodynamic tensor AI comes from solving for 
the Stokes flow around a particle moving with translational 
velocity U and and rotational velocity J7, assuming no slip 
on all particle and wall surfaces. 

If we then solve the electrokinetic problem for a fixed par- 
ticle in the applied field, we obtain the ICEO slip velocity 
Usiip as well as the total (hydrodynamic -I- electrostatic) force 
Fsiip and torque TsUp needed to hold the particle in place, 
thereby preventing ICEP and DEP motion. 

Using these calculations and invoking linearity, the condi- 
tion of zero total force and torque on the particle, 



'F slip 



0, determines the motion of the particle 



M 



slip 
slip 



(4) 



The particle trajectory is then described by the solution to 
the differential equation 



U 



together with the equations for the particle's angular orien- 
tation. 

This angular orientation does not matter for a full polar- 
izable or insulating particle. For the Janus particle, we will 
argue that only the rotations about x— axis are important, 
thus we will focus on the dynamics of just a single angle. In 
this case, the equation of motion is simply 



D. Nondimensional Equations 

We nondimcnsionalize the variables by 



X 



a EooO, 
eEr^a 



c 



EcqQ 



q = 



X 



D 

Ana 



L t" 



( V 



\eEl 



u =u 



\ V 



P 



Note that there are two time scales in the problem, r' = , 
the charging time, and t" = —Sr, the time scale for particle 
motion. 

Plugging in the equations, we obtain (after dropping the 
primes except for t) 

= 
V^u Vp 
Vu = 



with the boundary conditions 
dq 



dt' 

u = cv 



n ■ V(t) 



(5) 



on the particle surface, where C = 4>surface — 4'buik, is the zeta 
potential. For a polarizable particle, we have 4'surface = 
by symmetry, therefore we are left with C = —(j)buik — ~4'- 
In addition, we have 



V(/) X, |X| 



OO 



just as before, now with the dimensionless variables. The 
remaining equations are the no-slip boundary condition on 
the planar wall, and the requirement that the fiow vanishes 
at infinities. 

The constitutive relation between q and Q takes the simple 
form 

9 = -C = 

for the linear theory that we are going to analyze in this 
paper. 

The dimensionless force and torque on the particle are 
given by the formulae ^ and ([3]), where the stress tensors 
are replaced by their dimensionless counterparts 

GM = EE — 2^^-'- 

CTif = -pi + (Vu + (Vu)^) 



The force, angular momentum and stress tensors are scaled 
to 



?2 



Fref = sE^a, Tref = eE^a^ , aref = eE, 
Finally, the particle motion will be governed by 



dt" ' dt" 



E. Simplifications 

1. Steady Problems 

If a DC voltage is applied, then the system reaches a steady 
state after a while and the time derivatives drop out. This is 
the case when Neumann boundary conditions are valid also 
on the cylinder or sphere. In that case 



EE - ]^E'^l)ndA 



/( 

I {EE-n--E^n)dA = -- [ E'^ndA 
Jan 2 2 Jg^ 



because E • n =0 on the surface. As a consequence, the elec- 
trostatic torque induced on the particle is zero. 



2. Symmetry 

For the full cylinder problem, the electrostatic problem has 
an odd symmetry in x— direction, that is 

<i){x,z) = -<i}{-x,z) 

and for the full sphere problem, it has an odd symmetry in 
a;— direction and an even symmetry in y-direction: 

(t>{x, y, z) = (j){x, ±y, z) = -<t){-x, y, z) 

As a result, E^ has even symmetry in x and y. Therefore, 
in the steady case, the electrostatic forces vanish in those 
directions, and there can only be a vertical force. In general 
time-dependent cases, E has an even symmetry in x and and 
odd symmetry in y, and therefore the electrostatic force may 
not vanish in x direction. 

As for the Stokes problem, a glance at the slip formula 
shows that the slip, just like the potential field, has an odd 
symmetry in both x and y, and therefore so does the flow 
field. Consequently, there are no hydrodynamic forces in 
those directions. 

Needless to say, all these symmetry arguments disappear 
for the Janus particles. 



3. AC fields 

In the linear model wc arc considering, the time-periodic 
electrostatic forcing problem can be solved by letting 

(j) = Re(^e^'"*) 

and solving for the complex potential cf) using the equations 

V2^ = 
with the boundary conditions 

n ■ = iu}^ (polarizable) 
n • = (insulator) 

(f) = —EooX (at infinity) 

In the high frequency limit, the electrostatic problem ap- 
proaches to the solution of the dirichlct problem, that is, the 
first boundary condition is replaced by 



(polarizable) 



(6) 



This is because cj) — n ■ V4>/iuj ~ 0. Physically, this means 
that the double layers do not have enough time to charge 
when the forcing frequency is too high. 

Note that in cases of AC fields below, it is the field ampli- 
tude Eq that enters into the scalings above, and the U is the 
time-averaged velocity. 

Once the electrostatic potential is calculated, the time- 
averaged slip velocity can be obtained by the formula 



u, = -Re[C^;/] 



(7) 



where ( is the (complex) surface zeta potential, which is 
equal to —(j) in the linear theory, and E*^ is the complex 

conjugate of the tangential component oi E = W4>, the com- 
plex electric field. 

In the DC limit as w — > 0, the imaginary parts of the 

solutions go to zero, and we arc left with = ^C,E//^ which 
is the standard Smoluchowski's formula with a factor of 1/2. 



F. Numerical Methods 

We have solved the equations using the finite element soft- 
ware COMSOL for various geometries which have been ap- 
proximated by the cylindrical or spherical colloidal particle 
being enclosed by a finite rectangle. The equations are first 
converted to their weak forms, and entered into the general 
weak PDE module of COMSOL. None of the COMSOLs spe- 
cial modules are used. 
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FIG. 2: Streamlines to the Stokes flow problem for a cylinder 
near a wall. The surface color indicates the pressure. In this case, 
the increased pressure between the wall and the cylinder leads to 
repulsion away from the wall. 



For linear and nonlinear models alike, the computational 
efiiciency is improved by first solving the electrostatic prob- 
lem, and then the hydrodynamic problem. In time dependent 
cases, the fluid slip can be averaged and the Stokes problem 
is solved only once using this averaged slip. 

For reasons of completeness and easy reference, here we 
list the weak forms of the equations solved. This system 
of equations are converted to weak form by multiplying by 
corresponding test functions and integrating over the spatial 
domain. 

The electrical problem turns into 







n Jn 



In 



tqdr 



which is satisfied for all test functions (j>. The boundary con- 
dition for (f> is imposed in the form 

= y" <? C^coi - - C) rfr 

to be satisfied for all test functions q. 

The weak form for the stokes flow is similarly obtained as 

= - / [u • (V • a)dr+p\/ ■ u]dr 
Jn 

= [Vu : a—pV ■ \i\dv+ I u- (n • a) 
Jn Jan 



ds 



Since we do not have a simple expression for n • cr, it is best to 
introduce the new variable (Lagrangian multiplier) f = n • cr. 



This is also convenient for calculation of hydrodynamic forces 
at the surface. Then we obtain 



= - / [Vu : a-pV ■ u]dr~ 
Jn 



an 



[u • f + f • (u — Us)]ds 



III. ISOTROPIC PARTICLES NEAR A WALL 
A. Cylinder in a DC Field 

For isotropic particles near a wall, by symmetry, 
(^cylinder = 0, therefore = —(f). Moreover, there is no net 
horizontal force exerted on the particle, so the only forces of 
interest is in the vertical direction. Another consequence of 
symmetry is the absence of net torque on the cylinder. 

The DC cylinder problem has been solved analytically by 
Zhao and Bau [i^l in the linear case in bipolar coordinates. 
The mapping between the bipolar and the Cartesian coordi- 
nates is given by 



csin/? 



c sinh a 



coshf 



cos P 



, y 



cosh( 



cos/3 



where < a < oo, and — tt < /? < tt defines the region 
outside the cylinder. The geometric constants ao and c are 
defined as 



ao 



sech ^{a/h) 
h 



coth ao 



(note that there is a typo in the expression for ao in [47j). 
The hydrodynamic and electrostatic forces on the cylinder 
are calculated to be 



H 



27r sinh aoE'^c 



oo 

E 



(ao cosh ao — sinh ao) coth ao 
coshao 



1 



2 sinh ao 
1 



Fe 



^-^ V sinh(n -|- l)ao sinhao 

n— 1 ^ 



2TTEl.h 



coth ao 



sinh(rt 

n(n - 



2)ao sinhnao 
1) cosh ao 



}y 



sinh nao sinh nao sinh(n -I- l)ao 



Because of symmetry, there is no force in the horizontal di- 
rection. 

We are going to use this solution to gain some confidence 
in our numerical simulations. In Fig|31 you can see the com- 
parison of COMSOL results with the analytical expression. 
The match is especially good when the particle is close to the 
wall. It gets worse as this distance increases, because the ef- 
fects from the other walls also kick in. The simulation in 
Fig[3]has run with a box of size 20x20, and maximum mesh 
size 1, with finer mesh on the particle, specifically a maxi- 
mum size of 0.1. Experimentation with Comsol shows that 
the hydrodynamic error is sensitive to the size of the box, 
while the error in the electrostatic force is more sensitive to 
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FIG. 3: The Comsol numerical solution is compared to the ana- 
lytical one given by Zhao and Ban. Although the absolute errors 
tend to remain small, and the curves look identical, the relative 
error grows fast as the particle is located at larger distances from 
the planar wall. 



the mesh size. For a larger box, 40x40, and twice as finer 
mesh, we have obtained a similar picture with errors cut to 
about one third of their values in Fig[3] (results not shown). 



B. Cylinder in an AC field 

As the electric fields are screened quickly by the elec- 
trolyte, an AC field is usually preferred. Use of an AC elec- 
tric field also prevents harmful reactions on electrodes, and 
enables experimentalists to go to higher applied voltage dif- 
ferences. Such higher voltages may be desirable if they lead 
to stronger electrokinetic effects of interest. 

Far from the wall, the ICEO slip velocity around an ideally 
polarizable cylinder in an AC field was derived by Squires 
and Bazant [331, which takes the dimensionless form 



We) 



sin 26* 



(8) 



We use this expression to calibrate our numerical code, and 
find excellent agreement far from the wall. This result shows 
that ICEO fiow decays algebraically as above the RC 
charging frequency. Since electrostatic forces do not decay in 
this limit, we may expect a change in behavior near the wall. 
At high frequency, there is not enough time for double-layer 




total ^ ~ ' 



n 



150 



100 



50 



FIG. 5: Contour plot of total force on the ideally polarizable cylin- 
der. There is an equilibrium distance between the cylinder and 
the wall at high enough frequencies, indicated by the yellow con- 
tour line. As the frequency is increased, this distance decreases. 



relaxation, so the electric field ressembles that of a conductor 
in a uniform dielectric medium. 

An important observation is that the total hydrodynamic 
forces vanishes at higher frequencies whereas the total elec- 
trostatic force changes sign. As a result, if the frequency is 
high enough, there is an equilibrium distance from the wall. 
This distance decreases as the frequency is increased. 

In the high frequency limit, the electrostatic problem ap- 
proaches to the solution of the dirichlet problem, that is, the 
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laplace equation 



W^(j> = 



combined with the boundary condition 

(/) = (on the cyhnder) 
n ■ Wcj) — (on the wall) 

(f) — —EooX (at infinity) 

In this case, the solution is given by 

(j) = Re((^e'"*) = ^cosujt 

This problem can be solved analytically, and the solution is 
given by 



- 2cSoo 

OO 



OO r 



n=l 



cosh nao 

,-nao 



coshna sin nj3 



csin/3 



cosh nao 



cosh a — cos (3 
sin n/3 



Plugging this into the electrostatic force leads to the formula 



Fe,ui^oo = -2ttcEoo y^ ( 



J^::^ cosh nao 
n{n + 1) coshao 
sinh nao sinh (n + 1) ao 

with the same notation as in Zhao and Bau. 



Sphere in an AC field 



) 




FIG. 6: The hydrodynamic and electrostatic forces on a full metal 
sphere. 




(0/(X|^a/D) 



co/{X|^a/D) 



FIG. 7: Both the hydrodynamic and the electrostatic forces on 
an isotropic sphere are repulsive (away from the wall). While 
hydrodynamic forces on the sphere decline as a function of forcing 
frequency, electrostatic forces get stronger. 



ICEO flow around a sphere was first considered by 
Gamayunov et al. |1C| . Following the cylinder analysis of 
Squires and Bazant [39|, it is straightforward to derive the 
(dimensionless) ICEO slip velocity around an ideally polar- 
izable sphere in an AC field, far from the wall. 



{ue) 
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i2e 



16 1 + {uj/2f 



(9) 



Note that since (cos^ uot) — 1/2 the ICEO flow in a constant 
DC fleld Eq is twice as large as the time-averaged flow in an 
AC field E'o cos tot of the same amplitude: u^'~' = 2(u^'"}. 

For a sphere near a wall, the results are qualitatively the 
same as for a cylinder near a wall. That is, in the DC case, 
both the hydrodynamic and the electrostatic forces are re- 
pulsive. Moreover, the magnitude of hydrodynamic forces 
are about 2 orders of magnitude larger than the electrical 
forces. 

Note that the steady and time-periodic plots are consis- 
tent: At h/a = 1.5, the steady plots show forces about 
Fh ~ 8 and Fe ~ 0.12. The time-periodic plots, at the low 



frequency limit, start off at values Fh ~ 4 and Fe ~ 0.06, 
which are half of their steady counterparts. 

Unlike the cylinder problem, however, the electrostatic 
forces always remain repulsive, and therefore there is no equi- 
librium plane for the spherical particle. Instead, it is repelled 
to infinity by the wall regardless of the forcing frequency. 



IV. JANUS SPHERE NEAR A WALL 

A. Broken symmetries 

Without a nearby wall, a Janus sphere would align itself 
perpendicular to the electric field, that is, some of the elec- 
tric field lines would be included in the plane dividing the 
Janus particles metal and insulating sides. This effect has 
been demonstrated in (4Q] . and is presumably stronger than 
the wall effects, at least when the particle is sufficiently away 
from the wall. That being said, we will assume that the par- 
ticle always stays in the described configuration, that is, its 



dividing plane aligned with the electric field. This is not to 
say that the particle has no room for different rotational con- 
figurations, it can still rotate around x and y axis. Rotations 
about the y— axis (if existed) leave the particle unchanged, 
so we are left with rotations only around the x— axis. This is 
much easier to deal with than the original problem though, 
as only one angle is enough the describe the particles orien- 
tation. 

Far from the wall, the bulk velocity perpendicular to a 
DC field in the stable orientation is given by the formula of 
Squires and Bazant [l^l (Eq. 3.16), which takes the dimcn- 
sionless form, 



Udc 



9 

64 



2{Uac{uj ^ 0)) 



(10) 



neglecting compact-layer surface capacitance. As noted 
above, the time-averaged velocity in a sinusoidal AC field 
is smaller by a factor of two in the limit of zero frequency. 
Even in the bulk, without a wall, it is difficult to solve an- 
alytically for the ICEO flow at flnite AC frequency around 
a Janus particle, since the electrical response is not simply 
an induced dipole, due to the broken symmetry. Neverthe- 
less, we will see that the frequency dependence of the flow is 
similar to that around a sphere constant below the RC 
charging time and decaying above it. 

For a Janus sphere aligned perpendicular to the elec- 
tric field near a wall, a crucial observation is that the 
y— symmetry breaks down. As a result, there is a net force in 
the y— direction, as well as a net torque in a;— direction. The 
former leads to translation parallel to the wall, while the lat- 
ter causes rotation of the dielectric face toward the wall. We 
shall see that these effects of broken symmetry completely 
change the behavior near wall in an AC or DC field: Al- 
though a polarizable sphere is always repelled to infinity by 
an insulating wall, a Janus particle is always (eventually) 
attracted to it. 



B. 



Basic mechanism for wall attraction 



The key new effect is rotation due to hydrodynamic torque 
caused by asymmetric ICEO flow near the wall. This gen- 
erally causes the Janus particle to be attracted to the wall, 
as shown in flgure [5] The physical mechanism can be un- 
derstood as follows. When the fleld is first turned on, the 
Janus particle quickly rotates, by ICEP and DEP, to align 
its metal/insulator interface with the field axis, but with an 
arbitrary azimuthal angle, mainly set by the initial condi- 
tion. As described by Squires and Bazant [i^, the ICEO 
fiow around the particle draws in fiuid along the field axis 
and ejects it radially at the equator - but only on the polariz- 
able hemisphere, which acts like a "jet engine" drives ICEP 
motion in the direction of the non-polarizable hemisphere, 
which leads the way like a "nose" . 

Near a wall, as shown in the figure, the outward ICEO flow 
pushes down on the wall harder on the side of the polarizable 



u 




FIG. 8: Sketch of ICEO flows u and resulting ICEP torques T 
which cause a Janus particle to tilt its less polarizable end to- 
ward a wall, while translating toward the wall (until stopped by 
double-layer overlap) and perpendicular to the applied AC field 
E (directed into the page and parallel to the wall). This physi- 
cal mechanism may explain why the transverse ICEP motion of 
Janus particles was observable over the surface of a glass wall in 
the experiments of Gangwal et al. pd| ]. 



"engine" than on that of the non-polarizable "nose" , which 
produces a hydrodynamic torque tilting the nose toward the 
wall. A second cause of this rotation is the hydrodynamic 
coupling between ICEP translation parallel to the wall and 
rotation by shear stresses to cause rolling past the wall. Re- 
gardless of the initial position, these two sources of ICEP 
rotation cause the nose to eventually face the wall, so that 
the translational engine drives it toward the wall. This is 
likely the origin of the counter-intuitive attraction of Janus 
particles to a glass wall in the experiments of Gangwal et 
al [111]. 

What happens next depends on the details of the particle- 
wall interaction at very close distances. We will see that the 
bulk model with thin double layers must eventually break 
down, since the particle eventually collides with the wall, 
leading to overlapping particle and wall double layers. It is 
beyond the scope of this work to accurately treat the nonlin- 
ear and time-dependent behavior of these overlapping dou- 
ble layers, so we will explore two models: (i) infinitely thin 
double layers, i.e. using the bulk model to arbitrarily small 
heights, and (ii) a cutoff "collision" height, where overlap- 
ping double layers stop any further motion toward the wall, 
while still allowing transverse motion. The latter case as- 
sumes, as in the experiments TT], that the particles and walls 
have equilibrium surface charge of the same sign. For con- 
creteness, we will simulate Model (ii) with a cutoff height 
h = X = 0.05a, e.g. corresponding to a double-layer thick- 
ness (screening length) of A = 50nm with particles of size 
a = 1/im. 

Based on the simple examples above, we expect a subtle 
dependence on the AC frequency. Electrostatic DEP motion 
will always begin to dominate the hydrodynamic ICEP mo- 
tion at high frequency. Therefore, we now consider the low 
and high frequency cases separately. 
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FIG. 9: In the DC limit {uj ^ 0), we plot (a) horizontal velocity 
(b) vertical velocity and (c) tilting speed (degrees/charging time) 
as a function of the tilt angle 9 for the janus particle at distances 
h = 0.5a and h = 0.05a from the wall. 



FIG. 10: Typical trajectory of a janus particle under the DC limit 
uj —f interacting with the wall: As a function of time, plotted 
are (a) The horizontal speed (b) Distance from the wall (c) Tilt 
angle. Also, we plot the distance from the wall as a function of 
horizontal position in (d). 



C. Dynamics as a function of AC frequency 

As shown in Fig. ^ in the low frequency limit, the Janus 
particle experiences a rotational velocity turning its non- 
polarizable side toward the wall, as explained above. The hy- 
drodynamic ICEP torque is orders of magnitude larger than 
the electrostatic DEP torque, until the particle gets quite 
close to the wall. The magnitude of the horizontal ICEP 
velocity Uy parallel to the surface and perpendicular to the 
field is close to its bulk value Uy = 9/128 « 0.07 even fairly 
close to the wall at a height h = 0.5a at zero tilt, but reduces 
with the tilt angle. For small tilt angles and close to the wall 
a.t h — 0.05a, the horizontal velocity increases to Uy w 0.10, 
but it drops below the bulk value at larger tilt angles, e.g. 
to Uy « 0.05 at 6* = 45 degrees. Below we will see that this 
velocity is further reduced at higher forcing frequencies, due 
to the reduction of ICEO flow (since DEP cannot contribute 
to motion perpendicular to a uniform field). 

Regardless of the orientation, in the DC limit the particle 
moves ever closer to the wall in Model (i) since Uz < for 
any tilting of the nose toward the wall. Even if the the ver- 
tical motion is stopped at a critical height in Model (ii), the 
rotation continues in the DC limit until the particle points 
its non-polarizable nose directly at the wall {9 = 90) and the 
motion stops, although this can take a long time, since the 
rotation slows down substantially for tilt angles larger than 
45 degrees. As discussed below, a number of effects might 
lead to such a stabilization of the tilt angle, thus allowing 



steady translation along the wall. 

As shown in Fig. [TOl a typical simulated trajectory of the 
Janus particle shows it translating perpendicular to the field 
while rotating and attracting to the wall, until eventually 
coming to rest facing the wall. Even when the particle's 
motion stops, however, its polarizable hemisphere ("engine") 
continues driving a steady ICEO flow, which can lead to 
long-range hydrodynamic interactions with other particles. 
This is an interesting theoretical prediction which should be 
checked in experiments. Such immobilized Janus particles 
may have interesting applications in microfluidics. 

Similar behavior is predicted for finite AC frequencies in 
many cases. In particular, if a particle is initially mostly 
facing its non-polarizable hemisphere toward the wall {9 near 
90°), it will swim toward the wall and come to rest, as in the 
DC limit of Figure [TOl 

There are some new effects in AC fields, however, since 
ICEO flows are suppressed with increasing frequency. The 
competing effect of DEP can prevent the Janus particle from 
fully rotating and coming to rest on the surface, at least 
in Model (ii) where the collision is prevented artificially, as 
shown in Figure [T^l At w = 1 (the characteristic RC fre- 
quency of the particle), the rotation slows down substantially 
beyond 45° but does not appear to stop. In this regime the 
horizontal velocity decays to Uy ~ 0.015. For ut — 10 the 
particle appears to settles down to an equilibrium tilt angle 
around 45°, while steadily translating over the wall. The 
limiting horizontal velocity is roughly Uy ~ 0.009. 
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FIG. 11: For AC frequency ujtc = 1, we plot (a) horizontal ve- 
locity (b) vertical velocity and (c) tilting speed (degrees/charging 
time) as a function of the tilt angle 9 for the janus particle at 
distances h — 0.5a and h = 0.05a from the wall. 



FIG. 13: For AC frequency {ujtc — 10), we plot (a) horizontal ve- 
locity (b) vertical velocity and (c) tilting speed (degrees/charging 
time) as a function of the tilt angle 9 for the janus particle at 
distances h — 0.5a and h — 0.05a from the wall. 
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FIG. 12: Typical trajectory of a janus particle under AC fre- 
quency ujTc — 10 interacting with the wall: As a function of time, 
plotted are (a) The horizontal speed (b) Distance from the wall (c) 
Tilt angle. Also, we plot the distance from the wall as a function 
of horizontal position in (d) . 



D. Comparison to experiment 

The simulations with Model (ii) are in reasonable agree- 
ment with the experimental observations of Gangwal et 
al [ll| for metallo-dielectric Janus particles in dilute NaCl so- 
lutions in the low- frequency regime w < 1. The bulk theory 
of Squires and Bazant pTH) accurately fits the experimental 
velocity as a function of the field strength (Fig. 3 of Ref. [llj) 
and the particle size (Fig. 4), if a ratio 6 — 10 for the 
compact-layer to diffuse-layer capacitance is used to obtain 
the dimensionless formula, Uexpt = (9/128)/(l-|-10) = 0.006. 
However, the ICEP motion is observed only very close to the 
walls. 

Our simulations predict that the particles are quickly at- 
tracted to the walls over a time of order the channel width 
(60 /zm) divided by the typical ICEP velocity (10 /im/s), 
which is roughly one minute, consistent with experimental 
observations. The particles are also predicted to tilt, and 
moderate tilt angles can also be inferred from experimental 
images, although more accurate measurements are needed. If 
the tilt angle stabilizes around 45° (see below) , then the sim- 
ulations (Fig. [5]) predict that the ICEP translational velocity 
should be only 0.05/0.07 = 70% of the bulk value close to 
the wall, which would imply 6 = 7. This value is somewhat 
larger than that inferred from prior experiments on ICEO 
flow in dilute KCl around a larger (100/im radius) platinum 
cylinder j23|, but it is also observed that the ICEP velocity 
is slower than predicted at larger sizes (Fig. 4 of Ref. pH). 
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Apart from the rotational dynamics, therefore, the theory is 
able to predict the ICEP velocity fairly well. 

Without stopping the rotation artificially, we are able to 
predict the experimentally observed steady motion along the 
wall only at moderate to large uj. The reduction of ICEO 
flow in this regime reduces hydrodynamic torque (see below) 
and also enhances the effect of stabilizing electrostatic forces. 
Although Uexpt = 0.006 is measured in the low-frequency 
plateau < 1, this behavior otherwise seems quite consis- 
tent, since the slower ICEP velocity can also fit the exper- 
imental data using smaller (and perhaps more reasonable) 
values of 6. For example, the predicted velocity oiU = 0.015 
at (jj = 1 implies 5 = 1.5, while the velocity U = 0.009 at 

= 10 implies 5 — 0.5. 

The difficulty in predicting the stable tilt angle at low 
frequency may be due to our use of the low-voltage, dilute- 
solution theory, which generally overpredicts the magnitude 
of ICEO flows, especially with increasing salt concentration. 
For example, the electrophoretic mobility can saturate at 
large induced voltages, and the charging dynamics can also 
be altered significantly when crowding effects are taken into 
account 0. As a result, our simulation results at moder- 
ate frequencies lo — 1 — 10, which exhibit reduced ICEO 
flow due to incomplete double-layer charging, may ressemble 
the predictions of more accurate large- voltage, concentrated- 
solution theories at low frequency w < 1, where flow is re- 
duced instead by ion crowding in the double layer. This will 
be the subject of future work. 
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FIG. 14: Torque on a fixed Janus sphere versus tilt angle at 
heights h = 0.5o and 0.05a when (a) a; — > (b) ujtc = 10. 



E. 



Contact mechanics 




Another source of error in the model is our inaccurate 
treatment of the contact region, where double-layers overlap. 
We have simply imposed a small cutoff height in Model (ii) to 
prevent the wall collision, but there may be more complicated 
mechanical effects of the contact region. In particular, there 
may be enhanced hydrodynamic slip, due to the repulsion of 
overlapping (equilibrium) double layers of the same sign, as 
in the experiments. 

By examining the forces and torques close to the wall, we 
can infer to some degree what mechanical properties of the 
contact region might lead to the observed ICEP sliding along 
the wall and smaller tilt angles at lower frequencies (and thus 
also somewhat larger velocities). As shown in Fig.[T31 before 
the particle gets very close to the wall, the (mostly hydro- 
dynamic) torque acts to completely tilt the non-polarizable 
face toward the wall leading to collision. As noted above in 
Fig. [51 this can be understood as a result of the downward 
component of ICEO flow on the polarizablc hemisphere rais- 
ing the pressure by pushing on the wall on that side. 

The situation changes when the particle gets very close 
to the wall. As shown in Fig. [TH the torque changes sign 
at a tilt angle which is roughly 45°. This again can be un- 
derstood from Fig. [SJ since the ICEO flow between the par- 
ticle and the wall on the polarizable side, which drives the 



FIG. 15: Horizontal force on a fixed Janus sphere versus tilt angle 
at heights h = 0.5a and 0.05a when (a.) uj —> (b) lutc = 10. 



torque, is mostly absent. It would thus seem that even in a 
DC field, the particle would not rotate any farther, but this 
thinking neglects the hydrodynamic coupling between trans- 
lational force and rotational velocity near the wall, Eq. 
In Fig. [151 we see that the force on the particle parallel to the 
wall Fy remains strong, and this leads to a rolling effect over 
the wall due to shear stresses. For this reason, the rotational 
velocity persists in Fig. [HI even when the torque goes to zero 
in Fig. [ll 

The model assumes no slip on all non-polarizable surfaces, 
but this may not be a good approximation near the contact 
point when double layers overlap. If the equilibrium sur- 
face charges (or zeta potentials) on the non-polarizable hemi- 
sphere and the wall have opposite signs, then the overlapping 
double layers lead to a strong attraction, which would only 
stiffen the effective contact with the surface, and thus only 
increase the viscous rolling effect during motion along the 
surface. If the equilibrium surface charges (or zeta poten- 
tials) have the same sign, however, as in the experiments on 
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gold-coated latex Janus particles near glass walls [ll|, then 
there is a strong repulsion at the contact point. This repul- 
sion stops the collision with the wall in Model (ii) , but it may 
also "lubricate" the contact and allow for some sliding. This 
effective slip over the wall near the contact point could re- 
duce the viscous rolling, and, in the absence of torque, cause 
the rotation to stop, or at least be reduced for tile angles 
above 45°. In that case, we might expect a more accurate 
model of the contact region to predict to the experimentally 
observed motion, sliding over the surface by ICEP with a 
small tilt angle [9 < 45°), for a wider range of conditions, 
including lower AC frequency, perhaps even in the DC limit. 

V. CONCLUSION 

We have use the existing low-voltage theory of ICEP to 
predict the motion of polarizable particles near an insulating 
wall. Our results for symmetric spheres and cylinders con- 
firm the expected repulsion from the wall due to ICEO flow, 
sketched in Figure 1(a). In the case of the cylinder we show 
that attraction is also possible at high frequency, where DEP 
from electrostatic forces dominates slip-driven ICEP motion. 



Our results for asymmetric Janus particles reveal an un- 
expected attraction to the wall by a novel mechanism illus- 
trated in Figure [8l which involves tilting of the less po- 
larizable face toward the wall. Once it reaches the wall, if 
double-layer repulsion prevents further collision, the particle 
either rotates completely and ceases to move, while driving 
steady ICEO flow, or reaches an equilibrium tilt angle around 
45° while steadily translating along the surface, perpendic- 
ular to the elecric field. The latter motion only arises at 
moderate frequencies in our model, above the characteristic 
charging frequency for the double layers, while in experi- 
ments it is also observed at low frequencies. More accurate 
models taking into account reduced ICEO flow at large volt- 
age in non-dilute solutions and more accurate models of the 
contact region may improve the agreement with experiments. 

In any case, we have shown that polarizable particles can 
display complex interactions with walls due to broken sym- 
metries in ICEO flows. Attractive and repulsive interactions 
can be tuned by varying the geometry of the particles (and 
the walls), as well as the AC frequency and voltage. These 
phenomena may find applications in separations and self- 
assembly of colloids or in local flow generation in microfluidic 
devices. 
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